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Considering the collisional energy loss rates of heavy quarks from hard light parton interactions the 
total energy loss of a charm quark for a static medium has been computed. For the energy range E ~ 
(5 — 10) GeV of charm quark, it is found to be almost same order as that of radiative ones estimated 
to a first order opacity expansion. The collisional energy loss will become much more important 
for lower energy charm quarks and this feature could be very interesting for phenomenology of 
hadrons spectra. Using such collisional energy loss rates we estimate the momentum loss distribution 
employing a Fokker-Planck equation and the total energy loss of a charm quark for an expanding 
quark-gluon plasma under conditions resembling the RHIC energies. The fractional collisional energy 
loss is found to be suppressed by a factor of 5 as compared to static case and does not depend linearly 
on the system size. We also investigate the heavy to light hadrons D/tt ratio at moderately large 
(5 — 10) GeV /c transverse momenta and comment on its enhancement. 

PACS numbers: 12.38.Mh,24.85.+p 



I. INTRODUCTION 

In the initial stage of ultra-relativistic heavy-ion collisions energetic partons are produced from hard collisions 
between the partons of the nuclei. Receiving a large transverse momentum, these partons will propagate through the 
fireball which might consist of a quark-gluon phase for a transitional period of a few fm/c. These high-energy partons 
will manifest themselves as jets leaving the fireball. Owing to the interaction of the hard partons with the fireball 
these partons will lose energy. Hence jet quenching will result. The amount of quenching might depend on the state 
of matter of the fireball, i.e., quark-gluon plasma (QGP) or a hot hadron gas, respectively. Therefore jet quenching 
has been proposed as a possible signature for the QGP formation 1]. Indeed, first results from Au + Au at RHIC 
have shown a suppression of high_j3j_ hadron spectra spectra [2| which could possibly indicate the quenching of light 
quark and gluon jets [E 0> IE IS HI • O n the other hand the data from light ion interactions D + Au at RHIC 
indicate no evidence of suppression in high p± hadron spectra implying the absence of jet quenching as there is no 
formation of extended dense medium in the final state in such light ion interactions. However, this information from 
the light ion interactions in turn lends a strong circumstantial support that the observed suppression in Au + Au is 
due to the final state energy loss of jets in the dense QGP matter. 

Hadrons containing heavy quarks are important probes of strongly interacting matter produced in heavy ion colli- 
sions and has also excited a considerable interest. Heavy quark pairs are usually produced at early on a time scale 
of l/2Mc ~ 0.07fm/c from the initial fusion of partons (mostly from gg — ► cc, but also from qq — > cc) and also from 
quark-gluon plasma (QGP), if the initial temperature is high enough. There is no production at late times in the 
QGP and none in the hadronic matter. Thus, the total number of charm quarks get frozen very early in the history 
of collision which make them a good candidate for a probe of QGP, as one is then left with the task of determin- 
ing the p± distribution, whose details may reflect the developments in the plasma. The momenta distribution of c 
quarks are likely to be reflected in the corresponding quantities in D mesons as the c quarks should pick up a light 
quark, which are in great abundance and hadronize. The first PHENIX data |nj from RHIC in Au + Au collisions 
at y/s — 130AGeV on prompt single electron production are now available, which gives an opportunity to have an 
experimental estimate of the p± distribution of the heavy quarks. Within the admittedly large experimental error 
the data indicate the absence of a QCD medium effects. We hope that the future experimental study will provide 
data with improved statistics and wider p± range, which could then help us to understand the effect of medium 
modifications of the heavy quarks spectra. 

In order to see the effect of medium modifications on the final states, the energy loss of hard partons in the QGP 
has to be determined. There are two contributions to the energy loss of a parton in the QGP: one is caused by elastic 
collisions among the partons in the QGP and the other by radiation of the decelerated color charge, i.e., bremsstrahlung 
of gluons. The energy loss rates due to collisional scatterings among partons were estimated extensively 0, EH EI 
Ibi. Il5l He! IT7L Il8| in the literature. Using the Hard-Thermal-Loop (HTL) resummed perturbative QCD at finite 
temperature , the collisional energy loss of a heavy quark could be derived in a systematic way [2(J, lU |23, UM ■ 
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It was also shown 0, ^3 EB 13 that the drag force can be related to the elastic scatterings among partons in a 
formulation based on the Fokker Planck equation which is equivalent to the treatment of HTL approximations |2l| . 
From these results also an estimate for the collisional energy loss of energetic gluons and light quarks could be derived 
1251 . which was rederived later using the Leontovich relation 129, 

The energy loss due to multiple gluon radiation (bremsstrahlung) was estimated and shown to be the dominant 
process. For a review on the radiative energy loss see Ref. |2j|. Recently, it has also been shown 29] that for a 
moderate value of the parton energy there is a net reduction in the parton energy loss induced by multiple scattering 
due to a partial cancellation between stimulated emission and thermal absorption. This can cause a reduction of 
the light hadrons quenching factor as was first anticipated in Ref. Q, though the most of the earlier studies insisted 
that the collisional energy loss is insufficient to describe the medium modification of hadronic spectra. These studies, 
however, were limited to the case of massless energetic quarks and gluons. 

The first estimate of heavy quark radiative energy loss was found to dominate [3(| the average energy loss rate 
and subsequently the charmed hadron [3l| and the dilepton spectra had a strong dependence on the heavy quark 
radiative energy loss. Most recent studies of the medium modifications of the charm quark spectrum have computed 
by emphasizing only the energy loss of heavy quarks by gluon bremsstrahlung [33l IsiL Isfil l8fil ] . In Ref. it was 
shown that the appearance of kinematical dead cone effect due to the finite mass of the heavy quarks leads to a 
large reduction in radiative e nergy lo ss and affects significantly the estimation of the quenching of charm quarks and 
D/n ratio. Also in Refs. |34L l35l l36l| a surprising degree of reduction in radiative energy loss for heavy quarks was 
obtained by taking into account the opacity expansion with and without the Ter-Mikayelian (TM) effect. In the 
framework in which all these calculations have been performed the heavy quarks are possibly not ultra-relativistic 
(jv ~ 1) [l2l IrH fill l24l ] for much of the measured momentum range and in this case it is far from clear that radiative 
energy loss dominates over collisional case. 

Rest of the paper is organized as follows: in Sec. [H] the collisional energy loss for charm quarks is compared with 
the radiative ones computed in Refs. |34l l35j| . The quenching of hadron spectra in a medium is briefly reviewed in 
Sec. IIIII The charm quark in a thermally evolving plasma is modeled in Sec. HVI bv an expanding fireball created in 
relativistic heavy-ion collisions: we first obtain Fokkcr-Planck equation for a Brownian particle from a generic kinetic 
equation (Sec. IIV Af) : we, next, compute analytically the momentum loss distribution for charm quarks using the 
collisional energy loss rates through the transport coefficients based on elastic perturbative cross section implemented 
into a Fokker-Planck equation fSec. IIV B"|) ; the total energy loss for charm quark is obtained for an expanding plasma 
(Sec. IIV C|) : and then the quenching of hadron spectra and D/tt is obtained for RHIC energies (Sec. IIV D|) . We 
conclude in Sec. IVl with a brief discussion. 



II. HEAVY QUARKS IN A STATIC QUARK-GLUON PLASMA 



At leading order in strong coupling constant, a s , the energy loss of a heavy quark comes from elastic scattering 
from thermal light quarks and gluons. The energy loss rate of heavy quarks in the QGP due to elastic collisions was 
estimated in Ref. 21]. In the domain E « M 2 /T : it reads 



dE 
~dL 



ina x T 2 



1 



2v 2 



In 





In 









2^B(v) 

m 9 M 



(1) 



whereas for E » M 2 /T, it is 
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where nf is the number of quark flavors, a s is the strong coupling constant, m g — + rif /6)gT/3 is the thermal 
gluon mass, E is the energy and M is the mass of the heavy quarks. B(v) is a smooth velocity function, which can 
be taken approximately as 0.7. Following 11121) . one can now estimate the static energy loss for heavy quarks at the 
energies (temperatures) of interest. 

On the other hand, heavy quarks medium induced radiative energy loss |3^.l35| to all orders in opacity expansion, 
L/Xg (L is the length of the plasma, A g is the mean free path of the gluon), has been derived by generalizing the 
massless case to heavy quarks with mass in a QCD plasma with a gluon dispersion characterized by an asymptotic 
plasmon mass. This also provides the estimate of the influence of a plasma frequency cut-off on a gluon radiation 
(Ter-Mikayelian effect) and thus shields the collinear singularities (k± — > 0) those arise due to massless quarks. 

The medium induced radiative energy loss |35| for charm quark in first order opacity expansion, has been computed 
with a fixed Debye screening mass, /1 = 0.5 GeV with a s = 0.3, and a static plasma length, L = 4 fm with A g = 1 fm. 
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FIG. 1: Left panel: The scaled static energy loss of a charm quark AE/E as a function of energy E for a given length of the 
plasma, L = 4 fm. The collisional one is represented by solid line with plasma parameters for RHIC energy. The radiative 
energy losses according to Ref. [35| are also plotted in first order opacity expansion with (dotted) and without (dashed) the 
Ter-Mikayelian (TM) effect at a plasma length, L = 4 fm and a fixed Debye screening mass, fi = 0.5 GeV (see text for details). 
Right panel: The effective shift of the scaled collisional (solid curve) and radiative (dashed) energy loss AE/E as a function of 
distance L for a charm quark of energy E = 10 GeV. 



The scaled energy loss was found to obey a linear Bethe-Heitler like form, Ap:| lad « L ~ CL, where C is constant 

rad 

~ CE. Now, C can be estimated 



E 

d(AE) 



of proportionality per unit length. The differential energy loss follows as dL 

from right panel of Fig. ^ (also from Fig. 2 of Ref. [si|), as C ~ AJ ^" E ~ ^ip fm -1 at a plasma length L = 4 fm, 
/j, = 0.5 GeV and a s = 0.3. For a charm quark with energy E = 10 GeV, the differential radiative energy loss is 

estimated as ~ 0.375 GeV/fm. The Debye screening mass is given as /i = T Aira s {\ + ^f-) = 2.2415T, 

for two light flavors, rif = 2 and a s = 0.3. The Debye screening mass, fi — 0.5 GeV corresponds to a temperature, 
T = 0.225 GeV. With the plasma parameters corresponding to /i = 0.5 GeV, the differential collisional energy loss 

If A n\ COll 

for a 10 GeV charm quark in a static medium can also be estimated from 11121) as dL ~ 0.36 GeV/fm, and it 

is found to be of same order as that of radiative ones in Ref. [3{| . 

Now the total collisional energy loss can simply be evaluated from (11121) . The scaled collisional (solid line) and 
the radiative ones with (dotted) and without (dashed line) TM effect of a charm quark as a function energy E are 
displayed in the left panel of Fig. for a static plasma of length, L = 4 fm, with parameters T — 0.225 GeV (fi = 0.5 
GeV) and a s = 0.3. As discussed earlier that the radiative energy loss is proportional to E, resulting the scaled 
energy loss to be almost constant in E. In the energ y ra nge, E ~ (5 — 10) GeV, the scaled collisional energy loss is 
found to be almost similar as that of radiative ones [23 but decreases with E as the differential rates in (11121) have 
dependence on a log factor involving E. Therefore, the collisional energy loss will become much more important for 
lower energy range, and this feature itself will be quite interesting for phenomenology of particle spectra. 

In the right panel of Fig. we display the scaled effective energy loss of a charm quark due to collisional (solid 
curve) and radiative |3fij (dashed curve) ones in a static medium as a function of its thickness, L for a given charm 
quark energy E = 10 GeV. The thickness dependence of the scaled collisional energy loss for a given E is linear like 
the radiative case 0, whereas the earlier calculations [28l l37l l3S| show a quadratic form. This scaling clearly 
reflects a random walk in E and L as a fast parton moves in the medium 0, l35| with some interactions resulting in 
an energy gain and others in a loss of energy. 

In the energy range (5 — 10) GeV, which is much of the experimentally measured range (jyv < 4), the charm quark 
is not very ultra-relativistic and the collisions are found to be one of the the most dominant energy loss mechanisms. 
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In the weak coupling limit bremsstrahlung | 24| is the dominant energy loss mechanism if the charm quark is ultra- 
relativistic (p/v > 4). Though the collisions have a different spectrum than radiation, the collisional rather than 
radiative energy loss should in principle determine the medium modifications of the final state hadron spectra. In the 
following we would study the suppression of heavy quark spectra. 



III. QUENCHING OF HADRON SPECTRA 

We will follow the investigations by Baier et al. @ and Miiller Q , using the collisional instead of the radiative parton 
energy loss. Following Ref. 6] the p± distribution is given by the convolution of the transverse momentum distribution 
in elementary hadron-hadron collisions, evaluated at a shifted value p±^ + e, with the probability distribution, -D(e), 
in the energy e, lost by the partons to the medium by collisions, as 

da mcd f , , . da vac (p ± + e) f , , , dcr vac f , , . d da vac 
= / deD{e) -=ldeD(t)— + deD(e)e- 



d 2 pi_ J d 2 p± J d 2 p± J dp± d 2 p 



= d(jVaC | A£ d d<jVaC - da ™ c (P± + AE ) =Q( i3± ) dgVaC(:P±) (3) 
d 2 p± dp± d 2 p± d 2 p± d 2 p± 

Here Q(p±) is suppression factor due to the medium and the total energy loss by partons in the medium is 



AE = J eD{e)de. (4) 



We need to calculate the probability distribution, D(e), that a parton loses the energy, e, due to the elastic collisions 
in the medium. This requires the evolution of the energy distribution of a particle undergoing Brownian motion, 
which will be obtained in the following Sec. IIVI 



IV. CHARM QUARK IN AN EXPANDING PLASMA 

A. Generic Kinetic Equation, Fokker-Planck Equation, Drag and Diffusion Coefficients 

The operative equation for the Brownian motion of a test particle can be obtained from the Boltzmann equation, 
whose covariant form can be written as 

p^D(x,p,t) = C{D}, (5) 

where p^(E p , p) is the four momentum of the test particle, C{D} is the collision term and D(x, p,t) is the distribution 
due to the motion of the particle. If we assume a uniform plasma, the Boltzmann equation becomes 



3D _ C{D} _ (3D 



3t E V at 



(6) 



We intend to consider only the elastic collisions of the test parton with other partons in the background. The rate of 
collisions w(p, k) is given by 



(p,k)= ]T ^'(P,k), (7) 



w ^ 

3=Q,Q,9 

where w J represents the collision rate of a test parton i with other partons, j, in the plasma. The expression for uu J 
can be written as 

^(p,k)= 7j | ^^-(q)^, (8) 

where jj is the degeneracy factor, v TC i is the relative velocity between the test particle and other participating partons 
j from the background, Dj is the phase space density for the species j and cr? is the associated cross section. Due to 
this scattering the momentum of the test particle changes from p to p — k. Then the collision term on the right-hand 
side of © can be written as 



^ =J d 3 k [ W (p + k,k)^(p + k)- W (p,k) J D(p)] 



(9) 
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where the collision term has two contributions. The first one is the gain term where the transition rate w(p + k, k) 
represents the rate that a particle with momentum p + k loses momentum k due to the reaction with the medium. 
The second term represents the loss due to the scattering of a particle with momentum p. 

Now under the Landau approximation, i.e., most of the quark and gluon scattering is soft which implies that the 
function w(p, k) is sharply peaked at p rs k, one can expand the first term on the right-hand side of @ by a Taylor 
series as 

w(p + k, k)£>(p + k) « w(p, k)D(p) + k ■ ^-(wD) + hikj-^—iwD) + (10) 

op 2 opidpj 

Combining © , |J5J and (|10|) , one obtains a generic kinetic equation of the form 

^.|_|T„(p) D ] + s g-[ By( p) D ], (11) 
where the transport coefficients for momentum dispersion are given as 

T u (p) = Jd 3 kw(p,k)h = J d 3 kw(p,k) (p - p'), , (12) 

%(P) = ~ / d 3 kw(p,k)hk, = ~ / d 3 kw(p,k)(p~p') t (p-p'), . (13) 

These transport coefficients in H12I13|I depend on the distribution function, D, through the transition probability, 
w(p,k) in and can have different values depending upon the problem. The kinetic equation in (|lll) is the well 
known Landau equation |39j| . a non-linear integro-differential equation, which describes, in general, collision processes 
between two particles. It should therefore depend, in a generic sense, on the states of two participating particles in 
the collision process and hence on the product of two distribution functions making it non-linear in D. Therefore, it 
requires to be solved in a self-consistent way, which is indeed a non-trivial task. 

However, the problem can be simplified [39j if one considers a large amount of weakly coupled particles in thermal 
equilibrium at a temperature, T, constituting the heat bath in background and due to the fluctuation there can be 
some non-thermal but homogeneously distributed particles constituting foreground. It is assumed that the overall 
equilibrium of the bath will not be disturbed by the presence of such a few non-thermal particles. Because of their 
scarcity, one can also assume that these non-thermal particles will not interact among themselves but only with 
particles of the thermal bath in the background. This requires to replace the phase space distribution functions of 
the collision partners from heat bath appearing in © by time independent or thermal distribution, fj (q) . This will 
reduce the generic Landau kinetic equation l|llfl . a non- linear integro-differential equation, to Fokker-Planck (FP) 
equation, a linear differential equation for the Brownian motion of the non-thermal particles in the foreground. 

Now, one can write the transport coefficients in 1|12I13[1 for such a FP equation in terms of the two body matrix 
elements, M., between a foreground and a background particles [r2HT^ |: 

xbi-^]/(q)/(q) = (((p-p')i)) , (i4) 

#f/(p) = I{((p-P')i(p-P')j)) ■ (15) 

In our case, the incoming particle is a heavy quark which is different from the background. So, p(p') and q(q') 
represent the momenta of the incoming (outgoing) charm- and background light-quark/gluon, respectively. For each 
background species, there is a similar additive contribution to the collisional integral in l|14l) . 7 C is the spin and color 
degeneracy factor of the foreground particle arising due to initial reaction channels, /(q) is the particle distribution 
of the thermal background, and /(q) = [1 ± /(q)] corresponds a Bose enhancement/Pauli suppression factor for 
scattered background particles, as appropriate. Because of this thermal /(q) the contents of i|14ll5|) are different 
from (|12I13|1 and the charm quark in the foreground of a weakly coupled system is driven by a Brownian motion 
mechanism [H [H E HJ U E3 . 

We are now set to study the momentum distribution of a charm quark undergoing Brownian motion and its 
relation with the transport coefficients. In absence of vectors other than p the values of T^ p and Bj? . which depend 
functionally on p and the background temperature T, must be of the form in Langevin theory |l2Lll7l 39]: 

T* p (p,T) = Pl A( P 2 ,T) , (16) 
B?f(p,T) = Uj - B (p 2 ,T) + V -^B x {p\T) , (17) 
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where, p 2 = pf, summation convention is always implied. A is the drag, Bo is the transverse diffusion and ffi is the 
longitudinal diffusion coefficients. In terms of microscopic reaction amplitudes these functions are obtained [l2l Il7| 



A( P 2 ,T) = ((f)) - 



«p-p'» 



Bo( P 2 ,T) 



Bi( P 2 ,T) = - 



«(p-p') 2 » 



- 2((p-p')) +P 2 ((l)) 



(18) 
(19) 
(20) 



The averaging, ((• ■ ■ )), defined in Ijf 4|) can further be simplified |17| . by solving the kinematics in the centre-of-mass 
frame of the colliding particles, as 



«*V)» = 



f 



f 



51 2^ 4 7c E p 



\M\' 



q 2 rl 
-^q J-l 



d(cosx) 



V/( jBq ,)F(p') , 



4sM 2 



dcost 



(21) 



where Mc is the mass of a charm quark, s = (E p 



£ q ) 2 



(p + q) 2 , Eq' = E p + Eq — Ep/ and p' is a function of p, 



q and 9 c . m .. M 2 is the matrix elements [l2j for scattering processes Qg, Qq, and Qq, where Q is a heavy quark and 
17(g) is gluon(light quarks with 2-flavors). The expression in H21|) is larger by a factor of 2 than the ones originally 
derived in (3.6) of Ref. 1 2]. Apart from this we have also introduced the thermal masses of quarks (m g ), gluons (m g ) 
and the quantum statistics, as appropriate. 

The momentum and temperature dependence of the A, Bo and B\ are summarised in Figs. 1, 2 and 3 in Ref. fl?) ] 
and we do not repeat them here and refer the reader to this work for details. The main finding is that these coefficients 
are momentum independent upto p — 5 GeV/c and beyond this there is a weak momentum dependence. Note that 
the detailed studies of the dynamics of charm quark, as discussed in Refs. ma y om y depend on A and Bo, but 
perhaps not on B\ in the phcnomcnological relevant momentum range. In the present calculation we are interested 
in the kinematical domain of p — (5 — 10) GeV/c, for which one needs to solve the FP equation considering the 
momentum dependence of drag and diffusion coefficients. This will require to solve the FP equation numerically 

Instead, we assume the momentum independence of these coefficients in I|18I19I20|) . which will correspond to a 
scenario where a particle travels through an ideal heat bath and undergoes linear damping (Rayleigh's particle). So, 
these transport coefficients are expected to be largely determined by the properties of the heat bath and not so much 
by the nature of the test particle |39j |. This is also a fairly good approximation which we will justify in the next 
section. Under this approximation, the transport coefficients in I|16ll7f) become 



B 



FP 



FP 



PiA 
SuBo 



FP 



(22) 
(23) 



where B a (p -» 0,T) = 61 (p -> 0,T) = T 2 FP . This could be viewed as a course-grained picture in which the finer 
details of the collisions have been averaged out over a large number of macroscopic situation (or over an ensemble). 
Then combining (|llf> and (|22I23|) one can write the FP equation as 



at ~ fa [ u 1 



T-FP I 9 

2 w 



D 



(24) 



In the next Sec. II V lH we obtain the time evolution of FP equation in a thermally evolving QGP. 



B. Time Evolution of Fokker-Planck Equation, Drag and Diffusion Coefficients in an Expanding Plasma 

We assume that the background partonic system has achieved thermal equilibrium when the momenta of the 
background partons become locally isotropic. At the collider energies it has been estimated that to = 0.2 — 0.3 fm/c. 
Beyond this point, the further expansion is assumed to be described by Bjorken scaling law [43l | 

T(t)=tl /S T /tV 3 , (25) 
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FIG. 2: The momentum averaged ({A(t))) in 1311 and momentum dependence (A(p,t)) in 13UH of the drag coefficient of a 
charm quark in an expanding QGP with plasma parameters (see text) suitable for RHIC energy. 



where To is the initial temperature at which background has attained local thermal equilibrium. 
We consider, for simplicity, the one dimensional problem, for which FP equation in (|24|l reduces to 



dD d 



o 2 



dt dp 



(26) 



and as discussed in the previous Sec. II VAl that the coupling between the Brownian particle and the bath is weak, the 
quantities 7^ FP and 7^ FP can also be written using the Langevin formalism |39| as 



7f p (p) = J dkw(p,k)k= { -^ = (F)= P A 
1 ((Sp) 2 " 



T, 



2 St 

Now the work done by the drag force, 7^ FP , acting on a test particle is 

-dE = (F) ■ dL = T FP • dL , 
which can be related to the energy loss 0, |2{| of a particle as 

dE FP 



(27) 
(28) 

(29) 
(30) 



The drag coefficient is a very important quantity containing the dynamics of elastic collisions and it has a weak 
momentum dependence. Then one can average out the drag coefficient as 



(A( P ,T(t))) = A(T{t)) 



I dE 
p dL 



(31) 



implying that the dynamics is solely determined by the collisions in the heat bath and independent of the initial 
momentum of the Brownian particle. 

For averaging over the momentum the Boltzmann distribution and the differential energy loss rates I11I2H were used. 
The time dependence of the drag coefficient comes from assuming a temperature, T(t), decreasing with time as the 
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FIG. 3: Left panel: Comparison of diffusion coefficients given in 13 31 (solid line) and zero momentum limit of Bo (dashed line) 
given in If 91 . Right panel: Comparison of momentum averaged (T^ F ) in 134H and momentum dependence of Bo in If 9H . 



system expands, according to Bjorken scaling law given in 125(1 . We consider the initial temperature To = 0.5 
GeV, initial time to — 0.3 fm/c and a s — 0.3 of the plasma for RHIC energy. In Fig. [2] the momentum averaged as 
well the momentum dependence of drag coefficient of a charm quark in the QGP phase of the expanding fireball is 
shown as a function of the time. As can be seen that the behavior of the momentum averaged drag coefficient (solid 
curve) is dominated by T 2 /p ~ t -1 / 3 according to the scaling law. Now, it can also be seen that upto p = 10 GeV/c 
there is no significant difference between momentum averaged, (A(T(t))), and momentum dependence, A(p,T(t)), of 
drag coefficient and it has only a weak p dependence beyond p = 10 GeV/c. Since it decreases with moderately high 
values of p (> 15 GeV/c), the momentum averaged approximation of drag coefficient, (A(T(t))), would overestimate 
the actual A(p,T(t)) in this high momentum range. In our phenomenological approach the momentum independence 
of drag coefficient in ((31(1 is a good approximation upto a moderate value of momentum p < 15 GeV/c. 
Now, combining ((27(1 and 128(1 we can write the diffusion coefficient as 

T/ p =TT 1 FP =TA P . (32) 

Once the drag coefficient is averaged out using the properties of heat bath, one can approximate p by the temperature 
T of the bath (as discussed earlier that it is independent of initial momentum of the Brownian particle) and A by its 
average value given in ((31(1 . This leads 

T 2 FP =A(T(t))T 2 (t) , (33) 

which is also known as the Einstein relation [2l| between drag and diffusion coefficients. In the left panel of Fig. the 
diffusion coefficient obtained in ((33|) represented by solid line is displayed. It is found to have agreed quite well with 
the momentum independent diffusion coefficient (dashed line), Bo(p — > 0), in ((191) with a factor 1/3 multiplied with 
it, because we consider the 1-dimensional scenario. As evident that momentum independence of diffusion coefficient 
is also a fairly good approximation. 

Alternatively, one can also calculate the X, FP in 1(32(1 by substituting A from ((30(1 and averaging out the momentum 
dependence as 

<0 = r(-§). (34, 

In the right panel of Fig. [31 the average diffusion coefficient computed in ((34(1 (filled triangle) is displayed. The 
momentum dependence transverse diffusion coefficient, Bo(p 2 , T (t)), in I19|) is also plotted for different momenta. 
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It can be seen that there is weak momentum dependence in Bo in the momentum range p = 5 — 20 GeV/c. The 
momentum averaged values are in agreement with Bq for higher momenta p > 10 GeV/c whereas it overestimates 
lower momenta p < 10 GeV/c. We will use both the approximations for diffusion coefficient in jHSBH in our purpose 
to obtain the momentum distribution below. 
Combining lf2^1) and (|5Ujl . we find 



dD A d . ^ m d 2 D 



(35) 



where Dp, in general, has been used as the diffusion coefficient corresponding to (|33I34|) . and A is the averaged drag 
coefficient in 

Next we proceed with solving the above equation with the boundary condition 



D{p,t) 5(p-p ) . 

The solution of (|35|) can be found by making a Fourier transform of D(p, t), 



where the inverse transform is 



D(p,t) = — 
Zir 



D(x,t) 



+ 00 



D{x,t) e ipx dx 



D{p, t) e~ lpx dp 



Under the Fourier transform the corresponding initial condition follows from l|3tj[l and (|38[1 as 

D(x ,t = t ) = e - ipax ° 



(36) 



(37) 



(38) 



(39) 



where x = Xq at t = to is assumed. 



Replacing p 



i m and 



ix, the Fourier transform of (|35|l becomes 



9D t dD „ 2 ~ 



(40) 



This is a first order partial differential equation which may be solved by the method of characteristics [44j. The 
characteristic equation corresponding to l|40(l reads 



dt 

T 



Ox 
Ax 



dD 



V F x 2 D 



(41) 



Using the boundary condition in (|3l)|) the solution of l|4T)|) can be obtained as 



D(p,L) = 



VVW(L) 



exp 



P-Poe 



-SfA(t')dt' 



W{L) 



(42) 



where W(L) is given by 



W{L)=(aJ^ £> F (i')exp 2 J A(t")dt" dt'^j exp ^-2^ 



A(t') dt' 



(43) 



which is the probability distribution in momentum space. Since the plasma expands with the passage of time, we 
have used the length of the plasma, L as the maximum time limit for the relativistic case (jv ~ 1). 

In Fig. 0| we show the momentum loss probability distribution D(p, L) given in (|42|) . of a charm quark with initial 
momentum, po — 5 GeV/c, as a function of momentum p. The solid lines represent the distribution with the diffusion 
coefficient, T>p = AT 2 , in l|33|) whereas the dashed lines with (T>p) in l|34|) . Both set of curves are for two different 
expanded plasma lengths, L — 1 fm and 5 fm as indicated in Fig.0] In general, the physical mechanism reflected in 
this Fig. 01 can be understood at initial time (to) or length of the plasma a momentum distribution is sharply peaked 
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FIG. 4: The momentum-loss probability distribution, D(p, L) of a charm quark as a function of momentum p after plasma 
has expanded a distance L. 



at p = po, according to J3HJ). With passage of time (or distance traveled) the peak of the probability distribution is 
shifted towards smaller momentum, as a result of drag force is acting on the momentum of the charm quark indicating 
its most probable momentum loss due to elastic collisions in the medium. Moreover, the peak broadens slowly as a 
result of diffusion in momentum space, implying that a finite momentum dispersion sets in. As evident, with both the 
approximations in diffusion coefficient, only the momentum dispersion is affected while the peak positions remains 
unaltered, indicating that a drag force acting on the mean momentum of a charm quark is same. After plasma has 
expanded upto a length, L = 1 fm the charm quark loses 10% of its momentum whereas it is 25% at a expanded 
length of L = 5 fm. In the next subsec. llV CI we will use this distribution to compute the total energy loss of a charm 
quark for an expanding plasma. 



C. Energy-loss of a charm quark in an expanding plasma 

In the preceding subsec. TlV Bl we obtain a momentum-loss distribution by solving the time evolution of FP equation 
in a thermally evolving plasma, which is modeled by an expanding fireball under conditions resembling central Au-Au 
collisions at RHIC. The mean energy of a charm quark due to the elastic collisions in a expanding medium can be 
estimated as 

(E) = / ED(p,L)dp. (44) 
Jo 

The average energy loss due to elastic collisions in the medium is given by 

A£ = (e) = E Q - (E) , (45) 

where E = m± = \Jp\ + m 2 at the central rapidity region, y = 0. 

Using the momentum loss distribution in l)42|l the total energy loss of a charm quark has been computed in 145JI . 
The numerical results for scaled collisional energy loss of a charm quark as a function of energy E in an expanding 
plasma (solid line) is shown in the left panel of Fig.[5]for the plasma parameters Tq — 0.5 GeV, to = 0.3 fm/c, a s = 0.3 
and expanded plasma length upto L = 4 fm. In the energy range, E ~ 5 — 10 GeV, the fractional collisional energy 
loss remains almost constant around a value 0.15 and the reason for which can be traced back to the momentum 
independence of the drag coefficient 0> as discussed earlier. The corresponding scaled energy loss for a static 
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FIG. 5: Left panel: The fractional collisional energy loss of a charm quark AE/E as a function of energy E when plasma has 
expanded upto L — 4 fm (solid line) and for static plasma of length, L — 4 f m (dashed line). Right panel: Collisional AE/E 
as a function of length L for a charm quark of E — 10 GeV for expanding (solid line) and static plasma (dashed line) . 

plasma is shown by dashed line. Taking into account the expansion, the scaled energy loss is suppressed by a factor 
of 5 as compared to static case. In the right panel of Fig [5] the scaled energy loss is plotted as a function of length, 
L for a given energy, E — 10 GeV and found that it does not depend linearly on the system size for the expanding 
case (solid line) as compared to static case (dashed line). It is also expected that the similar suppression should also 
occur for radiative case due to the expansion |45|. 

D. Quenching of Hadron Spectra in an expanding plasma 

We assume that the geometry is described by a cylinder of radius R, as in the Boost invariant Bjorken model |43j 
of nuclear collisions, and the parton moves in the transverse plane in the local rest frame. Then a parton created at 
a point f with an angle <j> in the transverse direction will travel a distance 

L(<j>) = (R 2 - r 2 sin 2 </.) 1/2 - r cos0 , (46) 

where cos</> = v • f ; v is the velocity of the parton and r — \f\. The value of the transverse dimension is taken as 
R ~ 7 fm. 

The quenched spectrum convoluted with the transverse geometry of the partonic system can be written from @ 

as 

= Q(px) -a = TTTry) \ *t> / d r ^ " • ( 47 ) 



d 2 p±_ ^^ J <Pp±_ 2n 2 R 2 J r J d 2 p 

The p± distribution of charmed hadrons, D-mesons, produced in hadron collisions were experimentally found |4fj| 
to be well described by the following simple parameterization as 



d 2 Pl _ \bM 2 +p 2 ± 
where b = 1.4 ± 0.3, n = 10.0 ± 1.2 and M c = 1.5 GeV. 
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FIG. 6: Left panel: The ratio of charm to light quark quenching factors, Qh(p±)/Ql(p±), as a function of transverse 
momentum p± with collisional energy loss. Right pane: The ratio of charm to light quark quenching factors, Qh(p±)/Ql(p±), 
as a function of transverse momentum p± using both collisional and radiative energy losses. 



The parameterization of the p± distribution exists in the literature Q, IU E3 which describes the first RHIC light 
hadroproduction data for moderately large values of p± . In this case we consider the form given in Ref . [3J| which 
reads as 



dNl ac 



d 2 pi 



A 



1 



Po + P± 



(49) 



where m = 12.42 and po = 1.71 GeV/c. 

The light hadron quenching using collisional energy loss rate |25j was first anticipated || to be of same order as 
that of radiative Q ones, though most of the previous studies insisted that the collisional energy loss is insufficient to 
describe the medium modifications of hadron spectra. We also here estimate the light hadron quenching (for details 
see Ref. @) with the energy loss due to elastic collisions |2fj and the energy loss rate averaged over parton species 
reads 



dE 



+ + log 



2"// 2 ( 6 + ™/) Q Q2 . 



(50) 



We now illustrate in the left panel of Fig.[fj]the ratio of heavy to light quark quenching factors, Qh{p±)/Ql{p±), 
as a function of transverse momentum px, using the collisional energy loss for plasma parameters Tq — 0.5 GeV, 
to = 0.3 fm/c, a s = 0.3 and the charm quark mass, Mq = 1.5 GeV. As discussed earlier this ratio may reflect heavy 
to light hadrons, D/tt, ratio originating from the fragmentation of heavy and light quarks in heavy ion collisions. As 
shown the D/n ratio is enhanced significantly as compared to p — p collisions. The enhancement factor varies from 2.5 
to 4 in the p± range (5-10) GeV/c, due to the uncertainties of different choices of parameter set to parameterize the 
heavy hadron spectra as depicted in l|48|l . However, the ratio also strongly depends on the quenching of light quark 
jets. The numerical estimate shows that the quenching of charm quarks is about a half that of light quarks. The 
light quarks, for a given p±, lose 10% of their energy after traversing a distance 1 fm and around 40% after 5 fm Q 
whereas the charm quark loses 5% at 1 fm and around 20% at 5 fm (Fig. [SJ. Because of the large mass of the charm 
quark the D meson will be formed in shorter distance and hence charm quark would have less time to propagate in 
the medium before transforming into the D meson. On the other hand, the light quarks would travel in the medium 
over a longer period and suffer more loss in energy than heavy quarks. The ratio is also found to be little more than 
that obtained earlier |33| considering only the radiative energy loss with the appearance of the kinematical dead cone 
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effect due to the finite heavy quark mass. This implies that the collision is one of the most dominant mechanisms of 
energy loss in the medium. 

As shown that both collisional and radiative energy losses are of same order in magnitude it would be interesting to 
predict the D/tt ratio within our model considering radiative energy loss in addition to collisional one. Since neither 
the drag and diffusion coefficients have been calculated for other than collisional processes, it is not possible to infer 
the impact of radiative processes directly within our model. In order to circumvent our lack of knowledge of radiative 
processes in terms of the transport coefficients, we consider an alternative sets |48J within our model, obtained by 
multiplying the transport coefficients by a factor K = 2. It is our hope that the experimental data will allow us to fix 
an approximate value of K, if at all required. In the right panel of Fig.EJ we plot the D/tt ratio for such a case with 
the same parameter set as before. As shown the ratio is little reduced compared to the collisional case and it varies 
from 2 to 3 within the p± range, (5-10) GeV/c. There is no striking change in the ratio mostly due to the cancellation 
of the introduced K factor with reference to the collisional one for respective species. However, the small change can 
be attributed as the charm quark with the inclusion of radiative energy loss is relatively more quenched than that of 
the light quarks. 

One may also add an interesting scenario, after the quark jet has hadronized to leading particles, they would scatter 
with hadronic matter before decoupling. Considering udtt << 0W, e.g., it is most likely that the heavy mesons would 
decouple quickly from the hadronic phase. Pions would interact with each other via resonance formation and also 
with other light hadrons. This might lead to further enhancement of D/tt ratio. 

V. CONCLUSION 

Apart from uncertainties in the various parameters describing the plasma and hadron spectra let us also have a look 
at some of the assumptions made in this work which may affect our findings. First, as discussed above, the momentum 
dependence of the drag and diffusion coefficients, containing the dynamics of the elastic collisions, has been averaged 
out. A major advantage of this is the simplicity of the resulting differential equation. Of course, this simplification 
can lead to some amount of uncertainty in the quenching factor. Secondly, the entire discussion is based on the one 
dimensional Fokker-Planck equation and the Bjorken model of the nuclear collision, which may not be a very realistic 
description here but can provide a useful information of the problem. However, extension to three dimension is indeed 
an ambitious goal, which may cause that many of the considerations of the present work will have to be revised. 

In the present calculation the hadron spectra for both light and heavy quarks have directly been used to calculate 
the quenching factor, Q(p±). This is equivalent to assuming that a quark forms a hadron without much change in 
its energy. In order to calculate the effects of the parton energy loss on the quenching pattern of high p± partons 
in nuclear collisions, one should take into account the modification of the fragmentation function 0.0 of a leading 
quark resulting from many soft interactions of the hard partons in the medium. This also causes a significant energy 
loss of parton prior to hadronization and changes the kinematic variables of the fragmentation function [j| . This can 
also modify the quenching factor and thus D/tt ratio. 

We show that the total collisional energy loss is almost same order as that of radiative energy loss for a static 
plasma. Considering the collisional energy loss rates we obtain a momentum loss distribution for charm quarks by 
solving the time evolution of Fokker-Plank equation for a thermally evolving plasma. The total energy loss for an 
expanding plasma is found to be reduced by a factor of 5 as compared to static case and does not depend linearly on 
the system size. A ratio of heavy to light hadrons D/tt ratio is also estimated. Though the collisions have different 
spectrum than radiation and therefore contribute in a similar way to the suppression factor than anticipated earlier. 

We now eagerly wait for the experimental data and a detailed calculation have to be carried out before a realistic 
number for the value of the D/tt ratio can be presented. Nevertheless, the total collisional energy loss of a charm quark 
computed within our simplified model may imply that the collision is one of the important energy loss mechanisms 
in the medium for the energy range (5-10) GeV and this feature may be phenomenologically important. Also, the 
results for D/tt presented within this model is not definitive but can provide a very intuitive picture of the medium 
energy loss for partons in moderately large p± . 
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